DYNAMICS OF UNBINDING OF POLYMERS IN A RANDOM MEDIUM 
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We have studied the aging effect on the dynamics of unbinding of a double stranded directed 
polymer in a random medium. By using the Monte Carlo dynamics of a lattice model in two 
dimensions, for which disorder is known to be relevant, the unbinding dynamics is studied by allowing 
the bound polymer to relax in the random medium for a waiting time and then allowing the two 
strands to unbind. The subsequent dynamics is formulated in terms of the overlap of the two strands 
and also the overlap of each polymer with the configuration at the start of the unbinding process. 
The interrelations between the two and the nature of the dependence on the waiting time are studied. 
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I. INTRODUCTION 



The dynamics of polymers near a phase transition 
point, especially in presence of randomness or disorder, is 
important in many situations like denaturationjaf DNA, 
protein folding, collapse of heteropolymers etcEln. In ad- 
dition, polymers and extended elastic manifolds in ran- 
dom media constitute a class of problems which appear 
in various disguise in many problems, as e.g., interfaces 
in random systems, flux lines is high Tc superconduc- 
tors, surface growths and othersQ. Certain aspects of dy- 
namical properties for polymeric objects have been dis- 
cussed in the past, with emphasis mostly on the equi- 
librium or stationary dynamictu. In random systems, 
off-equilibrium dynamics has a special role because the 
system has to explore the phase space in search of its 
equilibrium state, if it reaches there at allQ. Thus, the off- 
equilibrium dynamics near a phase transition is expected 
to be different from the pure dynamics. In this paper, we 
study a very simple polymer model with a phase transi- 
tion for which equilibrium properties are known with a 
certain degree of confidence. The particular model we 
study is the unbinding transition of two interacting di- 
rected polymers in a random medium. This corresponds 
to a simplified model of denaturation af_DNA in a solu- 
tion with quenched random impuritiesB'El. 

Even though the off-equilibrium djynamics in glassy 
polymers are known for a long timdZI, the peculiarities 
of dynamics of random systems received attentioryjather 
recently through experiments on various systemsElu. As 
yet, there is no analytical approach for these problems, 
but several conflicting scenarios have been suggested, 
with the lack of well accepted equilibrium theories adding 
to the sore. In this respect, the model we are considering 
is in a rather enviable position, because of several analyt- 
ical tools and results available for equilibrium properties. 

A (d-l- 1) dimensional directed polymer (DP) is a poly- 
mer with a preferred direction so that it has random fluc- 
tuation in the transverse d directions only. Such an in- 



teracting DP system with homogeneous interactionJias 
been proposed in the past for denaturation of DNAQ in 
a pure solvent, where the most important degree of free- 
dom taken into account is the interstrand base pairing. 
Our model includes_a quenched distribution of impurities 
in the environmentEj. 

The maiiL effect of randomness in dynamics is the ag- 
ing effectd'til. If the system is allowed to equilibrate upto 
a certain time tw, to be called the waiting time, then 
the subsequent dynamics under a perturbation depends 
on this imposed time in a nontrivial way. We like to 
explore this aging effect in the dynamics of unbinding 
through a Monte Carlo dynamics of a lattice model. 

We discuss below in section 2 the equilibrium proper- 
ties of this interacting system of two polymers. There we 
also point out the connection of this two chain problem 
with that of a single chain via the replica approach. We 
then discuss the methodology of our simulation in section 
3. The results are presented and discussed in section 4. 



II. EQUILIBRIUM BEHAVIOR 

Let us consider two DP in the same random medium 
so that the Harniltonian in a path integral approach can 
be written asEJ 

i=2 „jv 1 

H = J2 dz[-Mzr + V{r,{z),z)] + 



N 



dz V ^(ri(2;) - r2(z)), 



(1) 



where r,i(z) denotes the d dimensional transverse spatial 
coordinate of the ith polymer at contour length z, and 
Vi{z) = dri{z)/dz. The first term denotes the elastic en- 
ergy part of the Gaussian chains and the second term 
is the random potential V{r,z) at point (r, z), and the 
last term denotes the mutual contact interaction between 
the chains. Note that the interaction is always at equal 
lengtlJa. 
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It is kipaai that randomness is relevant in d = 1 
dimensiont^H and the polymer has to swell to take ad- 
vantage of the favorable energy pockets. The transverse 
size grows with the length with an exponent = 2/3 
which is bigger than the Gaussian value 1/2 (expected 
for the pure case even in presence of the interaction. 

This particular problem of two interacting chain in 
the same random medium was considered numerically 
by Mezard in an attempt to calculate the overlap of two 
replicas for the single chain problem, the overlap ^ing 
the most important quantity in a replica approachtil. A 
general formulation for any d was given by Mukherji who, 
in addition to establishing the exact exponent for overlap 
in the 1-1-1 case, also obtained the relevant exponents for 
d = 2 + e for the spin glass transition point. This fpt- 
mulation was also used to study higher order overlapsEj, 
and in the strong coupling phasell3 for d > 2. 

In a .-dynamic renormalization group approach, 
Mukherjit2l showed that the interaction is relevant in all 
dimensions. Each chain individually behaves as in the 
single chain problem, i.e. the relevant strong disorder 
fixed point is independent oij-^j. A straight forward ex- 
tension of the approach of rellla gives the nontrivial fixed 
point for the two-jrepelling (i.e., unbound) chains in the 
random mediumllj. The fixed point diagram is shown in 
fig 1 , that shows that v — remains the critical point for 
the binding-unbinding transition for the two chains. A 
bound state forms for v < 0. The relevant exponents are 
also obtainable from the RG recursion relations. 

The order parameter that describes the critical point is 
the overlap or the number of contacts of the two chains, 
defined as 

1 

qiv) = j^ dzvSir,iz)-r2{z)). (2) 

The scaling behavior found for this overlap is g = 
f{vN2l^) for polymers of length N near the u* = fixed 
pointtHl. This particular scaling can be justified by a sim- 
ple argument. An overlap on a length scale L along the 
chain costs an energy vL while the gain from free energy 
fluctuation by following two different paths on this length 
L is L^, with X = 1/3 for this 1-1-1 dimensional problem. 
This gpes the scaling variable vl?'!^ as obtained exactly 
in Ref.EHl. We generalize this argument below for dynam- 
ics. The excitations we are considering here are the loops 
on a scale i|-and this forms the basis of the droplet pic- 
ture for DPEO. For large TV, the argument approaches 
the nontrivial fixed point, and being the unbound phase, 
(7 = 0, with a finite size scaling form q — /{AvN^^/"^). 
Av is the deviation from the fixed point. 

If we consider the single chain problem, then the over- 
lap, in the replica approach, is given by this q ai v — 0. 
Though this quantijUi, is not available from RG, numer- 
ical computationstSll3 show that g 7^ 0. This gives the 
Edwards- Anderson order parameter for this strong dis- 
order phase (see below). We therefore see that the order 
parameter for the critical point is a simple generalization 



of the order parameter needed for a replica approach of 
the single chain problem. In this respect, this DP prob- 
lem is unique among the known random models. 

In spite of these results for the equilibrium behavior, 
very little is known about the dynamics of unbinding, 
though certain aspects of, the. single chain dynamics have 
recently been looked intolijila. Our aim is to explore the 
time evolution of the overlap for the unbinding transition, 
and the effect of aging on this evolution, and correlate 
with the single chain behavior. 

III. MODEL AND METHOD 

To study the dynamics, we consider DPs on a square 
lattice. The polymers start at the origin and are allowed 
to take steps only in the +x or +y directions without 
any a priori bias. This produces polymers directed along 
the diagonal of the lattice. Two polymers interact if they 
share a point and each contact is assigned an energy u. 
In addition, there is a random energy at each site cho- 
sen from a uniform distribution e g[-.5,.5]. At a given 
temperature, there are two parameters,?; = u/ksT and 
e = e/ksT. We use the standard Metropohs single bead 
flip for the Monte Carlo dynamicsE2l. The chains^axe al- 
ways anchored at one end but free at the otheio. At 
each step the bead to be moved is chosen randomly from 
the 2N — 2 beads. One MC time step then consists of 
2N — 2 such attempts. The dynamics is performed for a 
given disorder realization, averaged over several random 
number realizations (thermal average)and initial config- 
uration, and then averaged over disorder realizations. 

Our procedure involves two chains completely bound 
(on top of each other) together evolving in the random 
potential for a time (i.e. MC is done with respect 
to random energy only) and then the chains evolve indi- 
vidually in presence of the interaction also. See Fig lb. 
With respect to the fixed point diagram of Fig la, the 
bound double stranded chain evolves towards the "strong 
disorder" fixed point K upto time t^, and after that the 
evolution is towards the stable fixed point M . We moni- 
tor the average fraction of contacts (overlaps) of the two 
chains and the overlap of each chain with the configura- 
tion at time i^,. 

Let us define two quantities self overlap Ci and mutual 
overlap q as 

C^{v,t + t^) = ^^{5{v,{t + U,)-v.,{t^))), 

q{v,t + t,,) = ^^{5{Ti{t + t^)-Y2{t + t^))), (3a) 

where (..) denotes thermal average and overbar denotes 
disorder average. The mutual overlap q{v,t) defined here 
is a time dependent generalization of the equilibrium 
overlap of Eq. H, while Ci is the overlap of the con- 
figuration of chain i at time t-\-tyj with its configuration 
at time tw By symmetry d is independent of the chain 
index i. 



2 



It is also possible to relate the overlap 9 to a corre- 
altion function. Let us define Si(t) = 1 if at time t 
there is an overlap of the two chains at chain length 
i, otherwise it is zero. The overlap at time t is then 
Si{t)/N . If we define an autocorrelation function 

C{ti,t2) = N^'^J2i{sii'ti)siit2)), we see that q{v,t + 
tw) = C{tw, t + tw) because Si{t^) = 1 forall i. 

For the single chain problem the self-over lap, Ci, de- 
fined above is also a quantity of fundamental importance. 
If we take limit t^, — > cxd first and then t 00, then for 
V — 0, Ci would correspond to the Edwards- Anderson 
order parameter for the strong disorder phase. This is 
because Ci would then measure the overlap, in equilib- 
rium, of the polymer configuration at two widely spaced 
time, and a nonzero value would imply a frozen random 
configuration, characteristic of a "strong disorder" phase. 
We therefore expect 

lim lim CiO,t + t^) = QEA- (4) 

t — ^00 tu^ — ^00 

In fact, for v — 0, in the limit i^, — > 00, one can also con- 
nect this overlap with the the self overlap defined above 
as q{0, t) = C(0, 2t) because in the equilibrium, the over- 
lap of the two configurations for C will be the same as the 
overlap needed for q. This is a check on our simulation 
for t ^tia. 

In the simulation, q and C were monitored for various 
values of v, and t^, for chains of length upto 300. At 
this length the dynamics we report here do not have sig- 
nificant finite size effects. Note also that by construction 
there is no finite size effect in the transverse direction. 



IV. RESULTS AND DISCUSSIONS 

We show the results of the simulation in Fig 2, where 
the overlap for various waiting times and v are plotted. 
Fig 2A shows the results for the pure system (e = 0) for 
V = 1.0, and there is no significant dependence on the 
waiting time. For the random case, shown in Fig 2B, we 
see the longer the waiting time the slower the relaxation. 
In other words, the system develops a stiffness as it ages 
in the random environment. This is the first effect of 
"aging". In absence of detailed theories, we considered 
various scaling forms. The form used for the single chain 
problem in Ref.E^ turns out to be applicable in this inter- 
acting problem. A data collapse is obtained by plotting 
cq(v, t + tw) vs. t/tw with suitable choices of the prefac- 
tor c. The variation of c with ty, and v is shown in the 
inset in Fig 2C. Similar scaling is obtained for the self 
overlap also (not shown). However, there seems to be no 
"universality" in the sense that the dynamics do depend 
upon the strength of interaction. It is not possible to go 
to large values of repulsion in this 1-1-1 dimensional prob- 
lem because of the log-jamming problem on a lattice. 

For t <^ t^, the early time dynamics is the "quasi" - 
equilibrium dynamics. For the largest t^^we find a linear 



relationship with the self overlap and the slope decreases 
with increasing v. Fig 3 A shows the early time equality, 
q{0, t) = C(0, 2t) for 1; = and its failure for v ^ 0. In 
fact, if we assume that for early times, t ^ t^,, C{t) « 
t-^, then, one can write q{0,t) = C{2t) = 2-^C(i). As- 
suming such a homogeneity relation for nonzero v, we 
can write q{v,t) = bC{v,t) so that by choosing the coef- 
ficient, b (= 2~^ for w = 0) it would be possible to get a 
data collapse for all v at least for early times. We do see 
such a collapse at early times as shown in Fig 3. This 
indicates a power law behavior, and we conclude that the 
early time power law decay of the overlap has the same 
exponent as the self overlap. 

Combining the various forms, a scaling formula for ag- 
ing can be suggested as 

q{v,t + t^,)^r''f{t/t^) (5) 

which for the limit — > 00, and then t — > 00, would give 
q — and not a finite qsA- Such a form has been used 
for various random systems in spite of this prohljinllil, 
and numerical simulations are yet to sort this onto. Fig 
3 suggests that a similar equation is valid for the self- 
over lapt^ with the same, rather small, exponent x. 

For the largest tn,, we see a power law decay of the over- 
lap at early times and not a logarithraic decay as would 
be expected from the droplet pictureEZl. In the droplet 
picture one assumes that the dynamics is governed by 
the typical barrier, and hence is of activated type. So, on 
a time scale i, the system would explore the phase space 
on length scales for which the barrier heights B Int. 
If one assumes further a growth of barriers with length 
scale B ^ L"^, then the relevant length scale at time t is 
L{t) ^ (Inty/^ . If we now generalize the scaling picture 
mentioned in section 2 to dynamics with the hypothesis 
that the dynamics is governed by the length scale L{t) at 
that time, one would expect a dynamic scaling 

q{v,t) = qivLitf'"^) = qiviXnti^-^^'"^) (6) 

This is valid for Int <C \ntw The simulation results are 
then not consistent with this dynamic scaling. In fact, 
no MC simulations have so far produced this log time 
scale in early dynamics in random systems. It has been 
speculated that the power law form, instead of logarithm 
of time, is a consequence of a logarithmic growth of bar- 
rier heights as opposed to B ~ . However, this is 
ruled out for Dlp-ibecause it is known from transfer ma- 
trix calculationsE2l that, in the 1-1-1 dimensional case, the 
typical barrier has the same scaling form as the free en- 
ergy fluctuation, = x = 1/3 . It is possible that the 
early dynamics is not controlled by the typical barriers 
but rather by the smallest barriers. If we denote the 
probability distribution of barrier heights by P{B), and 
if P{B) diverges (but integrable) as -B — > then early 
dynamics would not be activated type but rather like in 
spinodal decomposition where barrierless diffusion is the 
relevant mechanism. In such situations one finds j-that 
the time scale is a power law in the barrier hcighto as 
observed in simulations here. 
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In terms of lengthscales, the combined (bound) chains 
equihbrate by crossing barriers over lengthscales L{tyj), 
length being measured along the chain. The subsequent 
unbinding then involves the separation of the chains 
in presence of the repulsion within this length scale, 
L[t) <C Lityj). Once t !v t^,, one observes true nonequi- 
librium decay. Our data suggest again a power law but 
the overall decay of the overlap is rather small to get a 
reliable estimate of the exponent or any other functional 
form. However a scaling variable L{t)/L{tyj) seems to 
be a natural choice, which we find to be related to t/ty^. 
This also indicates that the relation between L{t) and t 
should be a power law type. We would like to add that 
there is the possibility that the length and time scales 
studied in these lattice simulations of this paper may not 
be in the right asymptotic limit to observe the dynamics 
predicted by the droplet picture. In fact, more analyt- 
ical work is necessary to understand the finite size and 
crossover effects in early dynamics of random systems in 
general. 

A bound on the late time decay of the overlap can 
be obtained by considering each bead independently (i.e. 
not connected as a polymer). In this case the overlap q 
is just the probability of reunion of two-vicious walkers 
(repulsive random walkers) at time nSa'Ea. This proba- 
bility for large times decays as , with = 3/2 for a 
pure system. Though its value for a random system is 
not known, it is expected to be smaller than the pure one 
due to the disorder induced effective attraction. For the 
polymer problem, the beads are connected and therefore 
this independent particle result gives an upper bound to 
the decay of overlap q for the polymer problem. The 
data for the pure system in Fig 2A can be fitted over the 
whole range by g = , with x = 1/3 < 3/2. 

The aging effects we have studied might be realized ex- 
perimentally also by letting DNA equilibrate in a random 
medium for a certain time and then suddenly changing 
the pH to start unbinding of the molecule. Early evo- 
lution of this unbinding will shed light not only on the 
dynamics of unbinding of DNA but also on the dynamics 
of random systems in general. 

In summary, we studied the aging effect of unbinding 
of a double stranded (directed) polymer where the fo- 
cus has been on the interchain interaction. The more 
time the double stranded molecule spends in the random 
medium the slower is the unbinding of the two strands. 
We have shown that the evolution of the overlap of the 
two chains has a scaling property where the time gets 
scaled by the waiting time in the random medium before 
unbinding. The average number of contacts of the two 
chains at early times evolve in the same manner as one of 
the strands as measured by the memory of its initial con- 
figuration, with a nonuniversal exponent that depends 
on the strength of the interaction. The late time decay 
that reflects the true nonequilibrium behavior shows also 
a power law behavior. Longer simulations are needed to 



clarify this nonequilibrium dynamics. 
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Fig 1 Bhattacharjee-Baumgartner 

FIG. 1. (a)Renormalization group fixed point and flow diagram for two cliains. G corresponds to tlie free pure gaussian 
polymers, F: pure, repulsive (fermionic or vicious walker) polymers, K: Strong disorder phase, M: repulsive polymers in random 
medium. The arrows indicate the flows of disorder and the interaction under renormalization.(b) The time sequence adopted 
in the simulation. 
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FIG. 2. (A) The overlap q(v, t + tn,) vs t for a pure system with w = 1. and f„ = 10, 10^, 10^, lO'*. (B) q{v, t + t-u,) vs t for 
V = 1.4 for a random system with = 10, 10^, 10^. (C) cq{v,t) vs i/i™ for (a) v = 0,(b) v = .2, (c) = .6, (d) v = 1.0, and 
(e) V = 1.4. The value of c is chosen for each data set (i.e., for each v and tw tw is taken as 10, 10^, 10^, 10* and 10^. The 
inset shows the variation of c with tw for t; = 0, .2, .6, 1., 1.4, v increasing upwards. 

o Fig 2 bhattacharjee-baumgartner 
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Fig 3 Bhattacharjee-Baumgartner 
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FIG. 3. Plot of q{v, t + tw) vs bC{v, t + tw) for t <tw, for various v and tw Inset A shows the plot of q{v, t + tw) vs C{v 
for the largest tw for each v (only a few data points are shown). The straight line is the equality line satisfied at w = 0. 



